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Abstract 



We report the spontaneous formation of a GaP fullerene cage in ab-initio 
Molecular Dynamics simulations starting from a bulk fragment. A systematic 
study of the geometric and electronic properties of neutral and ionized GaP 
clusters suggests the stability of hetero- fullerenes formed by a compound with 
zincblend bulk structure. We find that GaP fullerenes up to 28 atoms have 
high symmetry, closed electronic shells, large HOMO-LUMO energy gaps and 
do not dissociate when ionized. We compare our results for GaP with those 
obtained by other groups for the corresponding BN clusters. 

61.48.+C, 71.15.Pd, 81.05.Tp 



1 



Typeset using REVTpjX 



The discovery of carbon fullerenes and nanotubes |]||] has opened a completely new 
field at the borderline between chemistry and physics leading to many new phenomena and 
applications. 

Up to now, most efforts to identify fullerenes based on other elements have focused on 
BN which is the most similar to carbon and exists in nature in the hexagonal(graphite-like) 
structure [pHTOf] . However, the (nested-) cages and wires found for this material H|J do not 
resemble any of the small preferred structures of the carbon fullerene family, particularly 
due to the absence of the characteristic pentagonal rings. Besides, nanotubes based on other 
layered materials, such as GaSe [IT] and black-phosphorus [12], have been theoretically pre- 
dicted to be stable. On the basis of density functional calculations it has also been proposed 
that GaN nanotubes could be synthesised by using carbon nanotubes as a nucleation seed 

One intriguing question is whether fullerene cages could be realized in typical semicon- 
ductors of the III-V family, like GaAs, InSb or GaP, which do not possess a graphite-like 
bulk structure. These materials are not considered as good candidates for hollow structures 
since n bonding should be less effective in these larger atoms of higher rows of the periodic 
table than in the first one ||. 

In this letter we show, by means of ab-initio Car-Parrinello Molecular Dynamics |TJ]] , that 
a small GaP bulk fragment spontaneously organizes in a cage formed by a different number 
of atoms of the two elements arranged as in carbon fullerenes. We discuss the geometric 
and electronic structure of GaP cages with either the same or a different number of atoms 
of the two species. Our results strongly suggest that small GaP fullerenes could be stable, 
since they have high symmetry, closed electronic shells, large HOMO-LUMO energy gaps 
and do not dissociate when ionized. We give quantitative estimates of the relative stability 
of cages formed either by hexagons and pentagons as in carbon fullerenes or by hexagons 
and squares as proposed for BN [§-|9|] . 

Our results are obtained by the Car-Parrinello approach [jrjj using a Density Functional 
in the Generalized Gradient Approximation proposed by Becke and Perdew This 
approximation reproduces the experimental cohesive energy of typical bulk semiconductors 



within w 5% and underestimates the valence to conduction band excitation energies |T7 



We use nonlocal norm- conserving first-principles pseudopotentials JT8J and expand the single 
particle wavefunctions on a plane wave basis set with a cut-off of 12 Rydberg. We use a 
periodically repeated cubic simulation box of 24 A side, so that periodic images are at least 
14 A apart. We have verified that this size is large enough to describe isolated clusters. The 
electronic optimization and structural relaxation have been performed using damped second 
order dynamics with electronic mass preconditioning scheme [ITSL^D]]. We use throughout an 
integration time step of 8 a.u. The symmetry of the equilibrium structure is not biased but it 
is reached spontaneously during the geometry optimization starting from the corresponding 
regular polyhedron. 

The process of formation of the fullerene cage with 28 atoms from a larger bulk-like 
cluster of 41 atoms (Ga 28 P 13 ) is shown in Fig.l, with the help of three snapshots taken 
during the structural energy minimization which leads to the appearance of the Gai6Pi2 



fullerene cage |21|]. The cage has 12 pentagons and 6 hexagons and T^ symmetry as C28- An 
analysis of the charge distribution shows 108 valence electrons on the bonded cage, exactly 
the number which corresponds to the neutral Gai 6 Pi 2 cluster ]22J]. Finally in Fig. Id we show 
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the equilibrium structure of the neutral Gai6Pi2 cluster alone. 

The observed spontaneous formation of a Ga^P^ cage with pentagons is surprising 
since, in the case of BN [[)]-||, (deformed) squares are found to be energetically much more 
favorable. For B 12 N 12 ||, there is an energy difference of 9 eV between the cage with 
pentagons and the one with squares in favor of the latter which contains only heteropolar 
bonds and is favored for a material composed by atoms with very different electronegativity 
as B and N. Therefore, most studies have considered cages B n N n formed by hexagons closed 
by square rings ||||[8|||. Very recently, Fowler et al. |L(J have pointed out that, among the 
cages with pentagons, those with one species in excess of 4 atoms (B n N n+4 ) minimize the 
number of homopolar bonds. It is remarkable that the cage GaigPi2 which spontaneously 
appear in our simulation falls into this class. 

We have studied the equilibrium structure and electronic states of clusters with 20 and 28 
atoms of the type III n V n ± 4 , namely Gai 2 P 8 , Ga 8 Pi 2 , Gai 6 Pi 2 and Gai 2 Pi 6 , and compared 
them to clusters with the same number of III and V atoms, namely GaioPio and Gai 2 Pi 2 , 
the latter in the two isomers with hexagons closed either by five- or four-sided faces. 
The minimum energy structures of Gai6Pi2 and Gai 2 Pi6 are found to have T^ symmetry, 
whereas those of Gai 2 Ps and GasPi 2 have T^ symmetry. Among the clusters Ga n P n , the 
Gai 2 Pi 2 with 4-membered rings belongs to T^ whereas those with pentagons present very 
large distortions around the lower symmetry. The structural parameters of the cages 
belonging either to T^ or to T^ are given in Table I. Hexagons are found with alternating 
angles of 88° - 105° and 126° - 134°, pentagons with angles of 85° - 92°, 114° - 126°, ~ 100°, 
~ 110° and deformed squares with angles 75° and 98°. Ga-P distances are in general shorter 
than in bulk compounds due to predominant sp 2 bonding. The radial distance r from the 
center of the cluster given in Table I indicates a tendency of the anion to occupy positions 
at larger distances from the center than the cation, as found for ultrasmall cluster [^,^4]]. 

As in the case of carbon and BN fullerenes, the GaP clusters would represent metastable 
states with respect to the bulk equilibrium structure. Therefore only experimental obser- 
vation can establish with certainty their existence. Nevertheless, there are a few quantities 
which are used in the literature as indicators of stability. We support our prediction for the 
stability of the examined GaP clusters by using the following indicators: i) closed electronic 
shells and large energy gaps; ii) cohesive energy; iii) thermal stability; iv) stability of the 
ionized clusters. 

The first indicator of chemical stability is the energy gap between the highest occupied 
and lowest unoccupied molecular orbitals (HOMO and LUMO). In carbon fullerenes a cor- 
relation between this energy and the observed fullerenes has been experimentally verified 
| |25| |. In Table II we give the calculated HOMO-LUMO energy gap and the cohesive energy 
for all the clusters studied. Among the cages with pentagons the highest energy gaps are 
for cluster with P in excess of 4, a composition which has been suggested to be favorable 



also for BN [[nj. However, a very large gap is also found for the Ga i2 Pi 2 with squares. 

A comparison of the binding energies per atom between the GaP cages and the zincblend 
bulk phase of this material is possible only for the clusters with the same number of Ga and 
P atoms. From the results of Table II, we find that the cohesive energies per atom for 
Gai 2 Pi 2 with squares, Gai 2 Pi 2 with pentagons, and GaioPio are about 10% lower than in 
the bulk. This result is very close to that found for BN and carbon fullerenes of the same 
size |5|]. 
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We have studied the thermal stability of two clusters with very different energy gaps, 
namely Gai 2 Pg and GagPi 2 (see Table II). For both clusters we have performed two annealing 
cycles of about 3 ps, up to 1500 K and up to 2000 K. The system is heated with a rate 
of 2xl0 15 K/s, then equilibrated for one ps at the highest temperature, and finally cooled 
down with the same temperature change rate. For GagPi2 no bond breaking or structural 
rearrangements occur in both cycles and the structure comes back to the same minimum 
energy configuration when the temperature is lowered. This is also the case for Ga^Ps in 
the annealing up to 1500 K, whereas at 2000 K some structural rearrangement takes place 
leading to a distorted structure with higher energy when cooled down. These results indicate 
that the thermal stability is correlated with the width of the energy gap. 

Mass spectrometry experiments use the difference in mass-to-charge ratio of ionized 
atoms or clusters to select them. Therefore one basic requirement for the possible detection 
of such clusters is that they remain stable also when ionized. We have investigated the 
stability of some positively ionized clusters, [Ga 8 Pi 2 ] + , [Gai 2 Pi6] + and [Gai 6 Pi 2 ] + . We have 
included an uniform charge background in order to have an overall neutral system in the 
supercell calculation. The electronic structure remains almost unaffected and degeneracies 
are broken by negligible amounts in the order of hundredths of eV. Only minor structural 
distortions occur upon ionization. In particular, the six equivalent P-P (Ga-Ga) bond lengths 
split into three different classes. Remarkably, during a molecular dynamics run for [Ga 8 P 12 ] + 
we observe a dynamical exchange between these three classes of bond lengths with each other. 
This effect produces features in the low frequency vibrational spectrum in the range 30 — 120 
cm -1 which might be detected by infrared multiphoton ionization spectra . 



It is interesting to compare cages closed either by 4- or 5-membered rings. As already 
mentioned, such a comparison has been done for B 12 N 12 in Ref. @. However there are 
no results comparing clusters with equal number of atoms III n V n with square rings to the 
more favorable structures with pentagons of the type III n V n+ 4. The authors who have 
proposed the latter stoichiometry [|10| , in fact, do not give a comparison to cages with 
square rings. Although it is impossible to compare directly the cohesive energy of structures 
with a different number of atoms of each species, we are in a position to give an estimate of 
the bond energy for the two types of cages. Given the number of each type of bond in all 
the structures with pentagons studied so far and the values of the total cohesive energy, we 
estimate by a best fit the following energies per bond: Ec a ~p= -2.593 eV, Ec a -Ga=-1-133 eV, 



Ep_p=-2.349 eV ||27|| . As shown in Table II, these values yield the correct cohesive energy 
with a relative error of 0.4% at most for all cluster with pentagons, whereas overestimate 
the cohesive energy of the cluster Gai 2 Pi 2 with squares. In this cluster, in fact, there are 
only Ga-P bonds yielding directly Ec a _p=-2.499 eV a smaller value than in the clusters 
with pentagons. However, as it can be seen in Table II, it is the isomer with squares which 
has the lowest energy among the two Gai 2 Pi 2 . The 1.9 eV energy difference between them 
is much less than the 9 eV found for B 12 N 12 |J. This is most probably due to the less ionic 
character of the GaP bonds. The difference in electronegativity of Ga and P is in fact ~ 0.4 
against ~ 1 for BN. The observed spontaneous formation in our simulations of a cage with 
the same topology of C 2 s shows a possible evolution pattern from ionized bulk fragments to 
classical fullerene cages formed by pentagons and hexagons. 

In summary, we have shown, by means of ab-initio Car-Parrinello Molecular Dynamics, 
that small GaP fullerenes have highly symmetric structures, closed electronic shells and large 
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HOMO-LUMO gaps and cohesive energy. These clusters are thermally stable and remain in 
the same structure also when ionized. These findings together with the observed spontaneous 
formation in our simulations of a 28 atoms cage with the same symmetry of C 2 s support 
the possible existence of GaP fullerenes. We hope that our work will stimulate experimental 
groups to widen their search for hetero-fullerenes also to III-V compound semiconductors. 
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TABLES 
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Gai 2 Ps (Tft) 


Ga 8 P 12 (T ft ) 


Gai 2 Pi2 square (T h ) 


8 P 3 xxx x=2.05 r=3.55 

12 Ga i Oyz y=2.61 r=2.87 
z=1.19 


8 Ga 3 xxx x=1.54 r=2.67 

12 P i Oyz y=3.18 r=3.37 
z=1.13 


12 P i Oyz y=3.36 r=3.78 
z=1.74 

12 Ga i Oyz y=2.81 r=3.13 
z=1.39 



TABLE I. Point symmetry and atomic positions of the clusters with symmetry T^ and T^. 
Values of the parameters and of the radial distance r from the center of the cluster are in Angstrom. 
The positions of the minimum energy structure vary less than 0.01 A from the given values. 
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cluster 


HOMO-LUMO gap 

(eV) 


^cohesive 

(Hartree) 


E/it 
(Hartree) 


error 

(%) 


Gai 2 P 8 


1.28 


-2.5374 


-2.5370 


-0.01 


GagPi2 


2.14 


-2.8161 


-2.8052 


-0.39 


Gai6Pi2 


1.09 


-3.6876 


-3.6806 


-0.19 


Gai2Pi6 


1.55 


-3.9452 


-3.9488 


0.09 


GaioPio 


1.03 


-2.6632 


-2.6711 


0.30 


Gai 2 Pi2 (pentagon) 


1.24 


-3.2362 


-3.2429 


0.21 


Gai2Pi2 (square) 


1.86 


-3.3060 


-3.4308 


3.78 



TABLE II. HOMO-LUMO gap and cohesive energy of the clusters studied. For reference, the 
indirect energy gap of bulk GaP within the Local Density Approximation is 1.62 eV. The cohesive 
energy is obtained as the difference of the total energy minus the energy of the isolated pseudoatoms. 
Within the same approximations these values are -2.15045 Hartree and -6.46295 Hartree for Ga 
and P, respectively. The cohesive energy per GaP pair in the bulk is found to be 0.2976 Hartree. 
In the last two columns we give the values and the relative error of the cohesive energy estimated 
from the calculated bond energies Ec a -P= -2.593 eV, EG a -Ga=-L133 eV, Ep_p=-2.349 eV (see 
text). 
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FIGURES 



FIG. 1. Ga atoms are represented as large light balls, P atoms as small dark balls. Start- 
ing from a truncated bulk structure with tetrahedral symmetry, Fig. la shows a first step in the 
evolution towards structural energy minimization. The main rearrangement is the bonding of the 
peripheral Ga atoms between themselves, 12 atoms in pairs on the edges of the tetrahedron and 
four triplets on the vertices. In Fig. lb the central P atom breaks its bonds, followed by 12 Ga 
atoms (Figl.c) leading to the formation of a Gai6Pi2 fullerene cage. Fig.ld shows the equilibrium 
structure of the neutral Gai6Pi2 cluster alone. Notice the non-planarity of the pentagons. The 
symmetry of the equilibrium configuration is T^ (see Table I). 
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